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Abstract — Sun-tracking  microwave  radiometry  is  a  ground-based  technique  where  the  Sun  is  used  as  a 
beacon  source.  The  atmospheric  antenna  noise  temperature  is  measured  by  alternately  pointing 
toward-the-Sun  and  off-the-Sun  according  to  a  beam  switching  strategy.  By  properly  developing  an  ad 
hoc  processing  algorithm,  we  can  estimate  the  atmospheric  path  attenuation  in  all-weather  conditions. 
A  theoretical  framework  is  proposed  to  describe  the  Sun-tracking  radiometric  measurements  and  to 
evaluate  the  overall  error  budget.  Two  different  techniques,  based  respectively  on  elevation-scanning 
Langley  method  and  on  surface  meteorological  data  method,  are  proposed  and  compared  to  estimate 
the  clear-air  reference.  Application  to  available  Sun-tracking  radiometric  measurements  at  Ka,  V  and 
W  band  in  Rome  (NY,  USA)  is  shown  and  discussed  together  with  the  test  of  new  physically-based 
prediction  models  for  all-weather  path  attenuation  estimation  at  Ka,  V  and  W  band  from  multi¬ 
channel  microwave  radiometric  data.  Results  show  an  appealing  potential  of  this  overall  approach  in 
order  to  overcome  the  difficulties  to  perform  satellite-to-Earth  radiopropagation  experiments  in  the 
unexplored  millimeter-wave  and  submillimeter-wave  frequency  region. 


Index  Terms — Ground-based  microwave  radiometry,  Sun  tracking,  all-weather  path  attenuation, 
clouds  and  precipitation,  microwave  and  millimeter-wave  frequencies. 


I.  INTRODUCTION 

tmospheric  path  attenuation  is  one  of  the  most  important  limiting  factors  for  the  development  of  Earth- 
Asatellite  communications  at  Ka  band  and  beyond  [1],  [2].  At  these  frequencies  not  only  the  impact  of  rain  is 
significantly  affecting  the  channel  performances,  but  also  the  contribution  of  atmospheric  gases  and  non¬ 
precipitating  clouds  can  be  non  negligible  [3],  [4],  Radiopropagation  beacon  campaigns  at  Ka  band  and 
above  are  essential  to  experimentally  characterize  the  medium  behavior  at  these  frequency  bands  from  both  a 
physical  and  a  statistical  point  of  view  (e.g.,  [5]-[7]).  However,  these  campaigns  are  conditioned  by  the 
limited  number  of  dedicated  space  missions  or  by  the  opportunistic  satellite  constraints.  An  alternative  way 
to  approach  the  estimate  and  the  monitoring  of  path  attenuation  is  to  resort  to  ground-based  remote  sensing 
techniques  such  as  weather  radar  and  microwave  radiometers  (MWR)  (e.g.,  [8], [9]).  The  first  are  typically 
operated  at  frequencies  less  than  10  GHz  so  that  a  frequency  scaling  approach  is  needed  for  millimeter-wave 
estimates  [10].  The  latter  are  designed  up  183  GHz,  but  they  may  be  inadequate  in  presence  of  intense 
atmospheric  scattering  [11]. 

The  estimation  of  atmospheric  path  attenuation  from  ground-based  microwave  radiometers  (MWRs) 
may  be  affected  by  a  significant  underestimation  during  rain  events,  especially  at  frequencies  higher  than  Ku 
band  (e.g.,  [12], [13]).  These  errors  are  mostly  due  to  the  use  of  fixed  empirical  values  or  clear-air  based 
estimates  of  the  mean  radiative  temperature  within  path-attenuation  retrieval  algorithms  [14],  [15].  Model- 
based  approaches  have  been  proposed  in  the  recent  literature  to  overcome  this  problem  by  including 
multiple-scattering  effects  in  the  radiative  transfer  training  schemes  [16].  Indeed,  microwave  radiometric 
measurements  of  path  attenuation  during  rainy  and  cloudy  conditions  can  be  accomplished  in  a  quite 
accurate  way  by  exploiting  the  detection  of  the  solar  beacon  and  a  Sun-tracking  (ST)  operation  mode  [17]. 
This  approach  is  quite  well  known  in  radioastronomy  where  it  has  been  used  for  the  estimation  of  the  Sun 
brightness  temperature  at  microwave  and  millimeter-wave  frequencies  (e.g.,  [18]-[20]).  The  application  of 
ST  in  radiopropagation  for  estimating  attenuation  during  rain  was  envisaged  in  early  works  [21]-[23]. 

The  basic  idea  of  Sun-tracking  microwave  radiometry  (ST-MWR)  consists  in  using  the  Sun  as  a 
source  of  radiation:  the  sky  microwave  emission  is  measured  by  alternately  pointing  toward-the-Sun  and  off- 
the-Sun,  with  and  without  the  Sun  contribution  in  the  main  lobe,  according  to  a  beam  switching  strategy 
[17], [23].  Observations  toward-the-Sun  and  off-the-Sun  fulfill  two  main  goals:  i)  under  clear-sky  conditions 
they  allow  to  compute  the  Sun  brightness  temperature  at  the  different  frequencies;  ii)  during  rain  events,  they 
allowed  to  compute  the  path  attenuation  due  to  the  rainfall.  By  properly  choosing  the  switching  time  interval 
and  taking  into  account  the  main  lobe  aperture,  we  can  infer  the  atmospheric  attenuation  along  the  observed 
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path  in  all-weather  condition  through  indirect  evaluation  of  the  difference  between  the  two  measurements 

[17]- 

The  ST-MWR  technique  is  here  thoroughly  examined  and  applied  to  microwave  radiometric 
measurements  at  Ka,  V  and  W  band  [17], [24], [25].  The  available  dataset  consists  of  measurements,  collected 
in  2015,  by  a  ground-based  microwave  radiometer  with  a  tracking  control  stepping  motor  recently  installed 
in  Rome,  NY,  USA.  This  new  radiometer  has  four  channels  with  receivers  at  23.8,  31.4,  72.5  and  82.5  GHz 
and  is  an  adaptation  of  a  previous  system  [26].  The  exploitation  of  ST-MWR  in  this  experiment  is  aimed  at 
characterizing  the  millimeter-wave  atmospheric  channels  from  V  to  W  band  for  Geosynchronous  Earth-orbit 
(GEO)  and  Low  Earth-orbit  (LEO)  satellite  communications  without  resorting  to  ad  hoc  radiopropagation 
beacon  experiments. 

This  article  is  organized  as  follows.  In  section  II  the  theoretical  foundations  of  ST-MWR  will  be 
reviewed  and  main  assumptions  clearly  identified.  Subsections  II.A  and  II.B  will  discussed  two  different 
techniques  to  estimate  the  brightness  temperature  of  the  Sun  disk,  whereas  in  subsection  II. C  a  sensitivity 
analysis  framework  is  proposed  to  examine  the  predicted  performances  in  terms  of  errors  with  respect  to 
antenna  pattern,  beam  filling  and  Sun  brightness  temperature  uncertainties.  Section  III  is  devoted  to  the 
description  of  the  available  ST-MWR  data  and  to  the  development  of  physically-based  path-attenuation 
prediction  models  in  all-weather  conditions.  In  section  IV  the  comparison  of  ST-MWR  estimated 
atmospheric  path  attenuations  with  those  obtained  from  physically-based  parametric  prediction  models  at  Ka, 
V  and  W  band  is  presented.  Finally,  in  section  V  conclusions  and  future  developments  are  discussed. 

II.  Sun-tracking  microwave  radiometry 

For  ground-based  observations,  the  out-of-the-Sun  ( ooS)  sky  brightness  temperature  Tboos,  impinging  upon 
the  microwave  radiometer  antenna  along  the  zenith  angle  t?and  azimuth  (p,  is  given  by  [27]: 

TBOos(.0,V)  =  Tmrie,cp)[l  -  +  Tcose-^e^  (1) 

where  r  is  the  atmospheric  optical  thickness  or  path  attenuation  (in  Neper),  Tmr  is  the  sky  mean  radiative 
temperature,  and  Tcos  is  the  brightness  temperature  of  the  cosmic  background  (equal  to  about  2.73  K).  The 
frequency  dependence  of  parameters  is  here  neglected  in  favour  of  geometric  considerations.  The  antenna 
noise  temperature  Ta00s  along  the  antenna  pointing  angle  (S0,  <p0)  is  the  convolution  between  the  sky 
brightness  temperature  and  the  antenna  directivity  pattern  D{ Q,$)  expressed  by  [27]: 

Taoos(® O' Vo)  =  —  r^' B ooS  {d,(p)D{d0,(p0,e,cp)dn  (2) 

When  observing  the  Sun  with  an  antenna  solid  angle  Clant,  the  toward-the-sun  (twS)  antenna  noise 
temperature  Taiws  is  due  not  only  to  the  Sun  brightness  temperature,  but  also  to  the  sky  brightness 
temperature  emitted  by  the  observed  portion  within  the  antenna  beamwidth  so  that: 

TAtws(@o>(Po')  =  tz  In  TBtwS(d,(p)D(d0,(p0,d,(p')dn+—j  TBooS(d,(p)D(d0,(p0,d,(p)dn 

47T  “sun  47T  *Lsky 

(3) 

where  Qsty=  0.a„t-  Q™„  is  the  solid  angle  subtended  by  the  antenna  half-power  beam  subtracted  by  the 
Sun-disk  solid  angle  Clmn.  It  holds  that 


f„  D(9,(p)/DMdn  =  4tt 

XLant 


(4) 


where  Dm  is  maximum  directivity. 

The  twS  sky  brightness  temperature  is  given  by 
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TBtws(fi>  <P)  =  TBsune-T ^  +  Tmr{6,<p)[  1  -  e-«W]  +  Tcose~<9^ 

(5) 

where  7&„„  is  the  Sun  brightness  temperature.  If  we  can  approximately  consider  the  brightness 
temperatures  constant  within  the  beam,  then  (for  simplicity,  we  neglect  here  the  angle  dependence)  eq.  (3) 
can  be  rewritten  as 


Ta tws  ~  l’Btws  4n  D(0,  cp)d£l  +  TBooS  ^  fn,  D(8,  <p)d£l  —  TBtwS£l'sun  +  TBooS£ls 


471  n sky 


Lsky 


(6) 


where  the  prime  symbol  indicates  the  antenna  pattern  integration  effect.  If  £l'sky  =  £lant  ~  ^'sum 
substituting  the  brightness  temperature  expressions  and  dividing  by  flant,  we  get: 


Tam  =  fa[TBsune-T  +  Tmr(  1  -  e-)  +  Tcose~ T]  +  (1  -  /n)[Tmr(  1  -  e“T)  +  Tcose~ T] 


(7) 


being  the  beam-filling  factor  /n  the  ratio  between  the  directivity-weighted  Sun  solid  angle  and  the 
antenna  beamwidth  solid  angle  (see  also  (6)): 

r  _  flsun  _  ^-ant~^sky 

fn  n  n  ^  ' 

1 Lant  1 Lant 

By  reintroducing  the  angle  dependence,  the  previous  equation  for  the  twS  mode  can  be  further  rewritten  as: 
TAtwS  (0O. ^o)  =  faTBSune-r^^  +  Tmr(90,cp0)[l  -  C^o^o)]  +  Tcose~<9 


(9) 


Analogously,  for  the  ooS  mode,  under  the  same  assumptions,  it  holds 

TaoosV 0'<Pi)  =  TsoosiO 0,<Pi)  =  Tmr(60, <p1)[l  -  e~<9^]  +  Tcose~<9 ^ 


(10) 


The  different  azimuth  angles  <po  and  (pi  are  due  to  the  ST  operational  mode.  At  each  elevation  9o,  the  azimuth 
angle  is  switched  from  twS  <pi  to  ooS  (po. 

A.  Estimation  of  atmospheric  path  attenuation 

Assuming  that  the  switching  between  ooS  and  twS  observation  modes  is  fast  enough  and  the 
elevation  angle  9o  is  kept  constant  (so  that  atmospheric  attenuation  does  not  change),  the  Sun-tracking 
antenna  noise  temperature  difference  for  each  pointing  angle  can  be  expressed  by: 

ATa(90,  <Po’  <Pi)  —  TAtws(.@o>  <Po)  —  7aoo5(0o,  cpt)  (1  la) 

Note  that  the  beam-filling  factor  /n  can  be  dependent  from  the  pointing  angle  itself.  Substituting  (9)  and 
(19)  into  (11a)  and  assuming  the  mean  radiative  temperature  and  optical  thickness  do  not  change  between  the 
two  observation  modes  (that  is,  Tmr(90,  cp0)  =  Tmr  (0O,  <px)  and  t (0o,  cp0)  =  r(0o,  <Pi)),  we  get: 

ATA(90,(p0)  =  fn(90,(p0)TBsune-«9°^  = 

Tssune-^0^  (Hb) 
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From  (lib),  we  can  estimate  the  atmospheric  optical  thickness  x  (in  Np)  in  all-weather  condition  by 
means  of: 

=  <12a) 

and  in  terms  of  path  attenuation  A  (in  dB) 

A{G0,  cp0 )  =  4.343f  (0O,  cp0)  =  4.343 In  (12b) 

where  Tgsun  is  the  brightness  temperature  of  the  Sun  weighted  by  the  filling  factor  /n. 

Previous  expression  shows  that  path  attenuation  can  be  estimated  once  fa  and  weighted  Sun  brightness 
temperature  Tgsun  are  known.  The  first  one  can  be  derived  from  the  antenna  pattern  and  the  pointing  angle 
knowledge  (with  a  given  resolution  and  accuracy),  whereas  the  second  one  must  be  estimated  from  the 
available  measurements.  Note  that  the  dependence  of  the  beam-filling  factor  fa  on  the  pointing  angle  can  be 
neglected  if  the  radiometric  system  is  properly  designed  and  set  up;  an  evaluation  of  the  residual  errors  will 
be  shown  in  section  II. C. 

B.  Estimation  of  Sun  brightness  temperature 

As  mentioned,  in  order  to  estimate  atmospheric  path  attenuation  A,  the  weighted  Sun  brightness 
temperature  Tgsun  must  be  estimated.  Two  approaches  can  be  foreseen:  i)  the  Langley  elevation-based  self- 
consistent  method;  ii)  the  Tmr-based  meteorologically-oriented  method.  Both  methods  assume  the 
availability  of  radiometric  measurements  in  clear  air  conditions  where  a  plane -parallel  horizontally  stratified 
and  azimuthally  homogeneous  troposphere  can  be  assumed.  In  these  homogeneous  conditions  we  can 
hypothesize  that  Tgsun  estimates  are  minimally  affected  by  the  atmospheric  variability.  This  means  that  the 
“ secant  law ”  can  be  applied  to  describe  the  zenith  angle  dependence  of  antenna  noise  temperatures. 

(i)  The  Langley  technique  starts  from  (1 1)  which  can  be  rewritten  in  clear  air  as  (e.g.,[18],  [28]): 

ln[ATA(G0)]  =  ln[/nTBsun]  +Tc(r(0o)  =  ln[7^sun]  +  Tzc(r(0o)  sec(0o) 

(13) 

where,  due  to  the  previous  assumptions,  the  azimuth  dependence  has  been  removed,  and  the  following 
relation  holds  between  the  clear-air  slant  zclr  and  zenith  Tzclr  optical  thickness 

Tc;r(0  o)  =  Tzcir(0o)sec(0o)  =  Tzcir(0o,  <po)m(0o)  (14) 

In  (14)  m(0o)  stands  for  atmospheric  air  mass  and  is  equal  to  sec(0o).  By  plotting  ln[ATA(90)]  against 
air  mass  m(0o)  for  available  Sun-tracking  measurements  in  clear  air  and  approximating  the  measurements 
through  a  linear  best-fitting  curve: 


ln[ATA(90)]  =  a  +  b  m(0o)  (15) 

we  can  estimate  TBsun,  using  (lib),  from  the  intercept  of  the  previous  line: 

ffisun&o)  =  exp(a) 

1  Tzcir(0o)  =  -h  1  } 
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Once  Tgsun  is  estimated  during  clear-sky,  it  is  then  possible  to  retrieve  the  atmospheric  path  attenuation 
A(60,  <p0)  via  (12)  in  all  weather  conditions. 

(ii)  The  meteorological  technique  is  based  on  the  radiometer  equation  in  clear  air,  which  can  be 
written  from  (1)  as  [27]: 

TaciAB o)  =  TBclr(90)  =  Tmrclr[  1  -  e-rzcir(eo)/sin(0o)]  +  Tcose-^cir C«o)/sin(0o) 

(17) 


By  using  radiosounding  profile  measurements  of  meteorological  variables,  we  can  simulate  T'acIt  and 
Tmrcir ■  Indeed,  there  are  modeling  evidences  that  the  mean  radiating  temperature  Tmrdr  can  be  estimated 
directly  from  surface  meteorological  measurements  in  clear  air  (e.g.,  [29],  [30]): 

Tmrcir  ^0  d"  0-±Ts  +  Cl  2 Vs  ~f  Cl  3 (18) 


where  a,  are  the  regression  coefficients.  This  means  that  in  a  horizontally-stratified  clear  air  we  can 
obtain  the  zenith  atmospheric  optical  thickness  rzc;r(0o): 


Tzcir  (0o)  =  sin  (0o)/n 


Tmrcir  TCi 


1 


'-Tmrcir  TAcZr(^o)-l 


(19) 


From  (13),  using  the  Sun-tracking  measurements,  we  can  estimate  the  weighted  Sun  brightness 
temperature: 


ln[ATA(00)]  =  ln[/nTBsun]  +fcir(0o)  =  ln[7BsurJ  +  fzcir(0o)/sin(0o) 
that  is 


(20) 


fopflsun]  =  ln[ATA(60)]  -  fzcir(0o)/sin(0o) 

f  _  ffisun  (21) 

Bsun  fam 

Note  that,  with  respect  to  the  Langley  technique,  the  meteorological  technique  provides  a  time  series  of 
Tbsuu  on  a  clear-air  daily  basis.  The  retrieval  of  the  atmospheric  path  attenuation  A(90,  <p0)  via  (12)  is  then 
performed  by  taking  the  mean  value  of  fBsun  time  series. 

C.  Theoretical  error  sensitivity  analysis 

The  sensitivity  analysis  of  atmospheric  path  attenuation  estimate  to  residual  errors  or  uncertainties  of  ST- 
MWR  measurements  is  fundamental  to  understand  the  expected  accuracy  of  the  technique.  Several  sources 
of  uncertainty  can  be  identified  intrinsic  in  the  different  techniques  (Langley  vs  meteorological).  In  order  to 
perform  this  error  budget  analysis,  we  can  use  for  simplicity  the  first-order  error  propagation  theory  by 
assuming  a  statistical  independence  among  the  error  sources. 

Path  attenuation  is  directly  affected  by  antenna  noise  temperature  difference  uncertainties.  The  path 
attenuation  error  St  or  SA  (in  dB)  with  respect  to  the  antenna  noise  difference  uncertainty  SATa  is  obtained 
from  the  governing  equation  in  (12)  of  ST-MWR  leading  to: 


SA  =  4.3435r  =  -  ^  SATa  (22a) 

Fig.  la  shows  the  previous  equation  of  8 A  for  A TA  values  expected  in  clear-sky  situations  between  Ka  and 
W  band  (see  also  sect.  III).  As  ATa  is  typically  much  smaller  at  K  band  with  respect  to  V  band,  lower 
frequencies  are  more  affected  by  uncertainties  in  ATa  (due  to  calibration  errors,  antenna  mispointing  during 
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Sun  tracking  and  atmospheric  variability)  with  respect  to  higher  frequencies.  For  an  uncertainty  in  ATa  of  6 
K,  the  error  in  path  attenuation  goes  from  -0.05  dB  at  82.5  GHz  up  to  -0.27  dB  at  23.8  GHz. 

Fig.  lb  shows  the  same  analysis,  but  for  A TA  values  expected  in  rainy  situations.  Due  to  the  strong 
attenuation  of  Sun  emission  by  rain  layers,  ATa  differences  are  expected  to  be  much  smaller  than  those  in 
clear  sky.  For  heavy  rain  ATa  tends  to  zero  in  principle.  Although  the  sensitivity  increases  for  lower  ATa 
values,  possible  errors  associated  to  antenna  mispointing  of  the  Sun  position  are  less  important.  Nevertheless 
the  atmospheric  variability,  due  to  the  different  observation  geometry,  is  somehow  larger. 

As  noted  from  (12b),  since  ATa  values  may  be  negligible  for  heavy  rain  (when  the  Sun  disk  contribution  is 
completely  attenuated  by  the  atmosphere),  this  suggests  an  upper  limit  of  the  attainable  attenuation  during 
rain.  For  ATa  values  of  0.5  K  and  1  K  at  K  and  V-,W-band,  the  maximum  attainable  attenuation  with  the  ST 
method  is  about  25  to  28  dB,  respectively.  The  sensitivity  of  8A  with  respect  to  uncertainties  8TBsun  in  the 
estimate  of  TBsun  is  given  by 


8  A  =  4.3435t  =  ^  8T*sun  (22b) 

lBsun 

Fig.  lc  shows  8A  for  a  set  of  TBsun  values  which  are  those  expected  at  K  and  V  band  (see  sect.  III).  Alike 
Fig.  la,  lower  frequencies  are  associated  to  smaller  TBsun  values,  and  are  more  sensitive  to  errors  in  path 
attenuation.  For  an  uncertainty  of  TBsun  up  to  10  K,  the  error  in  path  attenuation  goes  from  0.06  up  to  0.36 
dB;  for  larger  uncertainties,  the  error  can  reach  1  dB  in  Fig.  lc. 


Fig.  1.  Sensitivity  analysis  of  ST-MWR  performances  for  a  set  ofTBsun  values  which  are  those  expected 
between  Ka  and  W  band.  Blue  line  corresponds  to  23.8  GHz,  red  line  to  31.4  GHz,  yellow  line  to  72.5  GHz 

and  finally  violet  line  to  82.5  GHz.  See  text  for  details. 

Note  that,  when  using  the  antenna  noise  temperature  difference,  we  have  assumed  the  same  elevation  both 
of-the-Sun  and  toward-the-Sun  observations.  If  this  is  not  the  case,  a  further  error  should  be  considered.  This 
means  that,  assuming  the  out-of-the-Sun  observation  in  (10)  is  performed  at  an  elevation  angle  9i,  the 
antenna  noise  temperature  difference  in  (1 1)  should  be  replaced  by: 


8Ta{Po,>  <Po<  ®i<  <Pi)  —  TAtws(.9o’  <Po)  TAoos(0  l’Vi)  (23) 


By  substituting  the  expressions  (9)  and  (10)  into  (23)  and  truncating  to  the  first  order  the  Taylor  expansion  of 
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the  atmospheric  transmittance  ratio,  we  can  obtain  the  following  TBsun  uncertainty  due  to  the  air  mass 
variation  8m  between  the  two  observations: 


STbsuu  =  ( Tcos  ~  Tmr)T2Clr  8m  (24) 

Fig.  Id  shows  the  previous  equation  for  a  set  of  tz  values,  which  are  those  expected  in  clear  air  between 
Ka  and  W  band,  and  for  a  fixed  mean  radiative  temperature  ( Tmr  =  270  K).  For  an  air-mass  uncertainty  of 
about  0.03,  corresponding  to  0.2  deg  at  20  deg  elevation,  the  error  in  estimating  TBsun  goes  from  -0.3  up  to  - 
2.5  K,  depending  on  the  value  of  the  considered  clear  air  optical  thickness  (higher  t zcir  values  give  higher 
errors). 


III.  Available  data  and  model-based  prediction 

This  section,  organized  in  2  subsections,  is  devoted  to  the  description  of  the  available  ST-MWR  data  and 
to  the  development  of  physically-based  prediction  models  at  Ka,  V  and  W  band  to  be  tested  by  using  ST- 
MWR  measurements. 

A.  Experimental  setup  and  MWR  data  processing 

The  available  dataset  consists  of  measurements  collected  by  the  ground-based  Air  Force  Research 
Laboratory  (AFRL)  MWR  during  May  2015  (20  days)  and  June  2015  (26  days)  in  Rome,  NY,  USA  (43.2°N, 
75.4°  W).  The  AFRL-MWR  has  four  channels  with  receivers  at  23.8,  31.4,  72.5  and  82.5  GFlz  and  is  a 
modified  version  of  the  RPG  LPW-U72-82  water-vapor  and  cloud-liquid  microwave  radiometer  series  in 
order  to  allow  an  automatic  Sun-switching  operation  mode  [26].  The  antenna  is  shaped  to  reduce  sidelobes  (- 
30dB  at  K-band  and  -40  dB  at  V-,W-band);  and  the  antenna  radiation  pattern  is  considered  approximately 
Gaussian  with  a  half-power  beamwidth  decreasing  with  the  frequency  and  equal  to  3.74°,  2.97°,  1.47°,  and 
1.30°,  respectively.  Measurements  have  been  collected  at  elevation  angle  between  20°  and  70°,  with  a  scan 
step  of  0.1°  in  both  elevation  and  azimuth,  during  clear  sky  days  and  during  raining  events. 

The  processing  and  quality-control  procedures  applied  to  the  AFRL-MWR  data  are  described  in  the 
following  paragraphs. 

1)  Clear-air  data  discrimination. 

For  selecting  clear-sky  days  in  order  to  estimate  TBsun  with  both  Langley  and  Meteorological  techniques, 
we  have  implemented  an  atmospheric  index,  a  scalar  quantity  named  Status  Sky  Indicator  (SSI),  purely 
based  on  the  available  radiometric  measurements  and  successfully  applied  in  several  MWR  campaigns 
[30] [31],  SSI  is  defined  as: 


SSI 


U?(31.4  GHz)  ~  c 
TB( 23.8  GHz') 


(25) 


where  c  is  a  parameter  dependent  on  air  mass  m.  If  SSI  is  minor  than  a  given  threshold  SSIth,  a  clear  air 
condition  is  assumed.  In  this  work  and  for  the  Rome,  NY,  site,  we  have  set  c  =  —0.13  ■  m2  +  6.3  ■  m  + 
2.1  and  the  SSIth  =  —0.00012  ■  m2  +  0.0066  ■  m  +  0.31. 

SSI  parameterization  has  been  set  up  by  performing  radiative  transfer  simulations  at  several  elevation 
angles  applied  to  a  long-term  available  radiosonde  observation  (RaOb)  dataset.  The  closest  RaOb  site  to 
Rome,  NY,  is  located  at  Albany  County  Airport,  NY,  USA  (WMO  station  ID  code  72518,  WBAN  ID  code 
14735).  RaOb  data  belonging  to  the  period  1994-2012  have  been  collected  for  this  study.  Downwelling 
brightness  temperatures  have  been  generated  using  a  plane  parallel  radiative  transfer  scheme  with  an  updated 
version  of  for  gas  absorption  and  cloud  model  [32]-[34], 
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2)  Filtering  toward-the-Sun  observations  in  clear  air  in  the  estimation  of  T*Bsun. 

In  order  to  estimate  TBsun  with  the  Langley  technique  (see  sect.  IV.A),  the  data  processing  has  been 
designed  to  select  the  maximum  TAtwS  for  each  elevation  angle.  The  AFRL-MWR  operational  mode 
maintains  a  constant  elevation  for  a  finite  time  so  that  the  Sun  can  move  inside  the  antenna  beam  during  that 
period.  As  noted  from  Fig.  la,  small  pointing  errors  of  the  Sun  position  can  lead  to  errors  in  A TA  of  several 
kelvins  with  a  large  impact  on  attenuation  estimates,  especially  at  K  band.  The  maximum  value  of  TAtwS 
ensures  to  identify  the  best  matching  observation  where  the  Sun  disk  is  centred  with  respect  to  the  antenna 
beamwidth. 

Finally,  a  binning  average  with  respect  to  air  mass  m  has  been  performed  considering  steps  of  0.2.  In  this 
way,  an  equal  distribution  of  samples  in  terms  of  air  mass  is  achieved  and  the  results  of  the  linear  regression 
in  (15)  are  not  influenced  by  the  different  distribution  of  samples  with  airmass. 


3)  Antenna  noise  temperature  self-consistency  check 

For  estimating  path  attenuation,  ST-MWR  approach  needs  valid  antenna  noise  temperature  differences 
A Ta  between  twS  and  ooS  modes  to  work  properly.  These  differences  may  reach  zero  or  even  negative  during 
intense  rain  events,  which  limits  the  application  of  this  technique,  as  detailed  in  Section  II.C. 

Only  differences,  greater  than  the  MWR  brightness  temperature  absolute  accuracy,  have  been  taken  into 
account.  The  latter  have  been  set  to  0.5  K  at  23-31  GHz  and  1  K  at  72-82  GHz  in  agreement  with  the 
manufacturer  specifications. 


23.8  GHz 


(a)  Samples 


(C)  Samples 


31.4  GHz 


(b)  Samples 

82.5  GHz 


(d)  Samples 


Fig.  2.  Time  series  of  ST-MWR  measurements  in  terms  of  ooS  (lower  couple  of  curves)  and  twS  (upper 
couple  of  curves)  antenna  noise  temperatures  for  two  case  studies  referring  to  a  clear  air  (8  May  2015,  blue 
dots  for  ooS  and  red  dots  for  twS)  and  cloudy  day  (28  May  2015,  green  dots  for  ooS  and  grey  dots  for  twS) 
and  for  the  AFRL-MWR  available  frequencies  a)  23.8  GHz;  b)  31.4  GHz;  c)  72.5  GHz;  d)  82.5  GHz. 


ST-MWR  measurements  are  interesting  to  analyze  with  respect  to  their  temporal  trend.  As  an  example 
here  discussed,  Fig.  2  shows  the  time  series  of  ST-MWR  measurements  in  terms  of  ooS  (lower  curves)  and 
twS  (upper  curves)  antenna  noise  temperatures  for  two  case  studies  referring  to  a  clear  air  (8  May  20 1 5)  and 
cloudy  day  (28  May  2015)  and  for  the  four  AFRL-MWR  available  frequencies  (Fig.  2a-d).  Fig.  3  shows  the 
same  as  in  Fig.  2,  but  for  two  different  case  studies  referring  to  a  moderate  rain  (1 1  May  2015)  and  intense 
rain  event  (30  June  2015). 

As  noted  in  Fig.  2a-b,  the  behavior  of  the  ooS  and  twS  antenna  temperatures  at  23.8  and  31.4  GHz  in  clear 
sky  is  similar  with  respect  to  the  elevation,  being  higher  values  observed  at  low  elevation  (at  the  beginning  of 
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the  daily  Sun-tracking)  and  reaching  minimum  values  at  the  solar  noon  (at  the  maximum  tracking  elevation). 
Conversely,  at  72.5  and  82.5  GHz,  the  behavior  of  TA  and  7^  time  series  is  the  opposite,  with  twS  time 
series  reaching  their  maximum  values  at  the  solar  noon  (Fig.  2c  and  Fig.  2d).  Such  behavior  is  explained  by 
recalling  (9)  and  (10)  and  the  different  impact  of  Tgsun  at  Ka  and  V,  W  band  (see  Fig.  lc).  At  Ka  band  the 
atmospheric  contribution  with  airmass  still  dominates  over  the  one  due  to  the  Sun,  whereas  a  V  and  W  band 
it  is  the  reverse.  In  presence  of  clouds,  the  antenna  noise  difference  between  the  two  measurement  modes 
ooS  and  twS  decreases  when  the  optical  thickness  increases  (i.e.,  for  each  frequency,  lower  curves  tend  to  to 
be  closer  to  the  upper  curves  in  Fig.  2).  This  is  especially  evident  in  the  rainy  event  where  the  difference  goes 
to  zero  as  atmospheric  extinction  significantly  increases  during  rainfall  (see  Fig.  lb  for  average  ATa  values 
during  rainy  conditions).  As  expected,  this  behavior  is  more  dominant  at  V  and  W  band  than  at  Ka  band. 


23.8  GHz 


31.4  GHz 


(a) 


Samples 

72.5  GHz 
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Fig.  3.  Time  series  of  ST-MWR  measurements  as  in  Fig.  2,  but  for  two  case  studies  referring  to  a 
moderate  rain  (11  May  2015,  blue  dots  for  ooS  and  red  dots  for  twS)  and  intense  rain  event  (30  June  2015, 
green  dots  for  ooS  and  grey  dots  for  twS)  for  the  four  AFRL-MWR  available  frequencies  a)-d). 

B.  Physically-based  prediction  models  of  path  attenuation 

The  ST-MWR  technique  can  offer  a  very  interesting  framework  to  validate  parametric  prediction  models 
at  frequency  bands  above  Ka  band  where  satellite-to-Earth  beacon  measurements  are  rare.  In  previous  works 
physically-based  prediction  models  (PPM)  have  been  proposed  for  estimating  specific  atmospheric 
parameters  (i.e.  attenuation  A,  Tmr  and  Tf)  as  a  function  of  selected  input  parameters  depending  on  both  the 
frequency  and  the  elevation  angle  [9],  [16].  These  models  have  been  based  on  the  non-linear  regression  fit  of 
numerical  simulations,  derived  from  the  sky-noise  Eddington  radiative-transfer  model  (SNEM)  in  an 
absorbing  and  scattering  medium  such  as  gaseous,  cloudy  and  rainy  atmosphere  [14],  [32], 

The  PPM  approach  is  based  on  the  exploitation  of  RaOb  datasets,  collected  in  a  location  close  to  the  site  of 
interest  and  used  to  statistically  characterize  the  local  meteorology  in  terms  of  temperature,  pressure  and 
humidity  average  and  standard- deviation  profiles  [24].  The  latter  statistics  is  then  imposed  in  the  Monte 
Carlo  pseudo-random  generation  of  vertical  cloud  structures  where  average  profiles  and  cross-correlation 
among  hydrometeor  concentration  are  imposed  [14].  The  considered  hydrometeor  categories  are  cloud 
droplets,  raindrops,  graupel  particles,  ice  crystals  and  snow  aggregates,  whereas  9  classes  of  cloud  structures 
(including  nimbostrati  and  cumulonimbi)  are  included  [16]. 

The  PPM  general  approach  has  been  adapted  for  Rome  (NY,  USA)  using  our  available  radiosonde  dataset 
from  the  site  of  Albany,  NY,  USA,  about  70  km  far  away.  By  extracting  the  meteorological  statistics,  SNEM 
simulations  have  been  performed  at  23.8,  31.4,  72.5  and  82.5  GHz  and  for  8  elevation  angles  between  20° 
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and  90°  in  order  to  compute  both  brightness  temperature  and  path  attenuation. 

Fig.  4  (a-d)  shows  the  scatterplot  of  TBooS  versus  corresponding  path  attenuation  A  at  zenith  (elevation 
angle  of  90°)  for  each  considered  frequency  in  all  weather  condition  (all  cloud  classes  plus  clear  air 
condition).  Fig.  5  (a-d)  shows  the  same  as  in  Fig.  4,  but  for  an  elevation  angle  of  36°.  Previous  scatterplots 
suggest  that  a  PPM  can  be  developed  to  estimate  path  attenuation  A  from  measured  antenna  noise 
temperature  (in  our  case  TBooS).  Two  parametric  non-linear  models  are  here  proposed  to  estimate  path 
attenuation  at  23.8,  31.4,  72.5  and  82.5  GHz  using  a  multiffequency  or  a  dual-frequency  MWR. 

The  first  PPM  model  is  a  multi-frequency  approach  and  it  is  based  on  a  polynomial  regression  (Pol)  on 
SNEM  dataset  [24],  [25].  The  proposed  multifrequency  PPM-Ro/  form  is  given  by: 

Apoiif)  =  m  ■  E?=  i  at  TBooS(fi )  +  bt  T$ooS(f{)  (26) 

where  /j=12, 3, 4=23. 8, 3 1.4, 72. 5, 82. 5  GHz  and  /  is  one  of  4  available  frequencies  f,  whereas  the 
coefficients  are  all  function  of  the  air  mass  m.  Table  A1  in  the  Appendix  provides  their  expressions  and 
numerical  values. 


m  tbN  id)  tb[ki 


Fig.  4.  Scatterplot  ofTBooS  versus  corresponding  path  attenuation  A  at  zenith  (elevation  angle  of  90°)  in 
all  weather  conditions  for  each  considered  frequency  a)  23.8  GHz;  b)  31.4  GHz;  c)  72.5  GHz;  d)  82.5  GHz. 


(0  Tb[K)  (d)  Tb[K] 


Fig.  5.  The  same  as  in  Fig.  4,  but  for  an  elevation  angle  of  36°. 

The  second  model  (PolDEx)  is  an  extended  version  of  PPM-Po/  approach,  reinforced  with  a  double 
exponential  single-frequency  term,  able  to  achieve  better  results  in  heavier  rainy  cases.  The  polynomial  and 
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the  exponential  terms  are  weighted  by  the  SSI  index  and  tuned  through  the  parameter  h.  The  PPM -PolDEx  is 
able,  substantially,  to  balance  heavy  weather  conditions  by  the  double  exponential  expression  with  the  good 
results  of  PPM-Ro/  in  cloudy  and  moderately  rainy  conditions.  The  multifrequency  PPM -PolDEx  expression 
is  expressed  by: 


ApoWExif)  =  MO-  -  SSI  +  K)APol (/)  +  ( SSI  -  h)ADEx(f )}  (27a) 

being  the  single-frequency  double-exponential  form  given  by 

ADEx(D  =  [ci 1  eC2'TBooS^  +  dx  ■  edrTBoosi.n ]  (27b) 

where,  as  in  (28),  the  coefficients  aL  and  bt  have  the  same  expressions  reported  in  Table  Al.  The 
coefficients  h,c1,c2,d1  and  d2  are  all  function  of  the  air  mass  in.  Table  A2  in  the  Appendix  provides  their 
expressions  and  numerical  values. 

To  highlight  the  advantages  of  our  all-weather  models,  we  can  also  show  the  comparison  with  the 
conventional  technique  using  the  clear-air  approximation  of  the  mean  radiative  temperature  Tmr  at  the 
frequency  / of  interest: 


AdAf)  =  4.343  In 


l 


Tmrclrif )  T cos 


■1 


L^mrcZrC/)  ^ AooSclri.^  o)J 


(28) 


where  Tmrcir  is  the  clear-air  mean  radiative  temperature  estimated  by  (18).  We  can  expect  a  significant 
error  when  trying  to  apply  (28),  instead  of  (26)  or  (27),  in  cloudy  and,  more  importantly,  rainy  conditions. 
Note  that  the  difference  Tmrcir  (/)  —  TAooScir(90)  is  less  than  zero,  the  estimator  in  (28)  is  not  applicable.  As 
mentioned,  we  have  conservatively  set  this  applicability  limit  not  to  zero,  but  to  0.5  K  for  23  and  3 1  GHz  to 
1  K  for  72.5  and  82.5  GHz. 

In  order  to  evaluate  the  accuracy  of  Pol  and  PolDEx  PPM  prediction  models,  we  can  evaluate  the  error  of 
the  regressive  parametric  estimates  with  respect  to  corresponding  SNEM  simulations.  Table  A3  in  the 
Appendix  shows  the  intercomparison  between  Pol  and  PolDEx  models  in  terms  of  root-mean- square-error 
(RMSE)  and  average  error  (AvE),  where  the  error  is  defined  as  the  difference  between  the  considered  model 
and  the  simulation.  We  can  clearly  note  from  Table  III  how  the  PolDEx  parametric  model  shows  better 
results  at  all  frequencies  and  for  all  elevation  angles.  RMSE  and  AvE  values  increase  with  the  decrease  of 
the  elevation  angle,  and  this  behavior  is  expected  as  optical  thickness  (and  path  attenuation)  increases  with 
the  decrease  of  the  elevation  angle  thus  meaning  that  the  percentage  error  is  still  comparable  at  any  angle. 

In  order  to  stress  the  latter  consideration,  an  index  of  agreement  (IA)  has  been  also  considered  to  better 
evaluate  the  percentage  accuracy  [35].  IA  is  a  standardized  measure  of  the  degree  of  model  prediction  error 
and  it  varies  between  0  and  1.  An  agreement  index  score  of  0  suggests  no  agreement  between  the  PPM  model 
and  the  SNEM  dataset,  while  an  agreement  score  of  1  suggests  complete  match  between  the  model  and  the 
dataset.  IA  is  defined  as  [35]: 


IA  —  ^ _ <  ( XPPM~XSNEM )2) _  ( 29 \ 

<  (  \xppm~{xsnem)\  +  frsNEM-(*SNEM)l 

where  xPPM  and  xSNEM  are  the  PPM  estimate  and  SNEM  reference  simulation.  Considering  the  PolDEx 
model  at  zenith  (elevation  angle  of  90°)  in  Table  III,  we  can  note  as  IA  is  equal  to  about  0.99  at  23.8  GHz 
and  about  0.84  at  82.5  GHz.  The  accuracy  decreases  with  the  decrease  of  elevation  angle  but  the  IA  is  always 
greater  than  0.7,  that  is  PPM  -PolDEx  is  able  to  provide  reliable  results  with  respect  to  SNEM  dataset. 

Worldwide  distribution  of  MWR  with  channels  in  V  and  W  bands  is  still  very  limited.  For  this 
reason  a  further  dual-frequency  polynomial  model  has  been  developed  in  order  to  provide  a  stand-alone 
solution  only  using  channels  at  Ka  band,  i.e.  at  23.8  and  31.4  GHz  (e.g.,  [6]).  Indeed,  in  order  to  retrieve  path 
attenuation  estimate  for  a  dual-channel  MWR  we  might  use  the  double-exponential  DEx  model  in  (27b) 
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since  it  is  single-frequency  and  SSI  exclusively  depends  on  the  brightness  temperatures  at  23.8  and  31.4 
GHz.  However,  an  extension  to  a  dual-frequency  frequency  can  improve  the  results.  The  dual-frequency 
PPM -Pol2  model  is  given  by: 


ApoM  =  m  ■  Z?=i  ai  TB00S(fi )  +  bt  TiooS{fd  (30) 

where  /i=i,2=  23.8,31.4  GHz  and  the  coefficients  are  all  function  of  the  air  mass  m.  Table  A4  in  the 
Appendix  provides  their  expressions  and  numerical  values.  Properly  replacing  APoi(f )  in  (27)  by  ZlPoi2(/) 
given  by  (30),  the  PolDEx  model  can  be  used  to  retrieve  path  attenuation  in  all-weather  conditions  from  the 
MWR  measurements  at  23.8  and  31.4  GHz  only. 


IV.  Application  to  sun-tracking  data 

This  section  is  devoted  to  the  description  of  results  obtained  from  ST-MWR  in  terms  of  path  attenuation 
estimates  at  Ka,  V  and  W  band.  This  work  has  been  performed  according  to  the  following  steps:  i)  careful 
selection  of  clear-sky  days  with  minimum  atmospheric  variability;  ii)  computation  of  TBsun  by  means  of  both 
Langley  and  meteorological  methods;  iii)  application  of  the  Sun-tracking  technique  to  estimate  total 
atmospheric  attenuation  in  all  weather  conditions  with  a  focus  on  rainy  cases;  iv)  application  of  PPM  Pol  and 
PolDEx  prediction  methods  to  estimate  all-weather  total  attenuation  from  ooS  observations;  v)  inter¬ 
comparison  analysis  of  the  retrieved  attenuations  from  steps  iii)  and  iv).  Verification  of  PPM  Pol  and 
PolDEx  prediction  algorithms  is  also  discussed  together  with  an  error  budget  analysis. 

A.  Results  for  Sun  brightness  temperature  estimate 

By  using  the  antenna  noise  temperature  difference  versus  air  mass,  Fig.  6  (a-d)  shows  the  estimate  of 
TBsun  using  the  Langley  technique,  as  discussed  in  sect.  II.B,  for  each  frequency  and  two  selected  clear-air 
days  (May  8  and  May  21,  2015).  The  best-fitting  line  is  also  shown. 

Data  processing  steps  for  clear-sky,  as  described  in  III.A,  have  been  performed  to  obtain  robust  estimates 
of  Tpsun  m  order  to  ensure  as  much  as  possible  a  constant  daily  optical  thickness.  For  the  Langley  technique, 
this  means  to  obtain  a  linear  dependence  with  a  relatively  small  deviation  in  (15),  being  rzcir  the  slope  in  the 
linear  approximation. 

TABLE  I.  Mean  radiative  temperature  regression 

COEFFICIENTS  FOR  APPLYING  THE  METEOROLOGICAL  TECHNIQUE 


1 

/[GHz]  | 

«ollG  | 

%  [K/K] 

■ 

<*2  1 
[K/mbarl 

a3  [K/%] 

23.8 

-118.65  | 

0.962 

0.107  ! 

0.146 

31.4  ! 

-125.65  ! 

0.99 

0.102  ! 

0.162 

72.5 

-110.88 

0.932 

0.102 

0.155 

82.5 

-134.49  | 

1.03 

0.099  | 

0.170 

By  using  the  times  series  of  antenna  noise  temperature  difference,  Fig.  7  (a-d)  shows  the  estimate  of  TBsun 
using  the  meteorological  technique  for  each  frequency  and  the  two  selected  clear-air  days.  The  average  value 
of  TBsun  is  also  highlighted  as  a  constant  line.  The  regression  coefficients  ao,  ai,  a2  and  a,?,  used  to  estimate 
Tmrdr  in  (18),  are  reported  in  Table  I  for  completeness. 

In  terms  of  Sun  emission  variation  at  the  involved  frequencies,  we  have  assumed  that  the  long-term 
component  is  small,  also  considering  that  the  antenna  beamwidth  is  much  larger  than  the  Sun  arch.  Indeed, 
we  have  not  observed  noticeable  variations  of  the  estimated  Sun  brightness  temperature  in  our 
measurements,  even  though  this  should  be  proved  on  a  longer  dataset. 
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Fig.  6.  Estimate  ofT^sun  using  the  Langley  technique ,  as  discussed  in  sect.  II.B,for  each  frequency  and 

the  2  selected  clear-air  days. 
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Fig.  7.  Estimate  ofTgsun  using  the  meteorological  technique  for  the  two  selected  clear-sky  days  (May  8, 
blue  dots;  May  21  cyan  dots)  at  each  frequency  a)  23.8  GHz;  b)  31.4  GHz;  c)  72.5  GHz;  d)  82.5  GHz. 


Table  II  shows  the  intercomparison  between  the  Langley  and  meteorological  technique  estimates  for 
T(isun  together  with  the  confidence  intervals  (in  terms  of  standard  deviation  for  the  meteorological  technique 
and  the  maximum  deviation  for  the  Langley  technique).  Note  that  Table  I  also  provides  the  values  of  t zcir 
and  the  maximum  deviations  oTzclr  (which  is  below  5%  at  all  frequencies  thus  proving  a  good  selection  of 
clear-sky  days).  The  uncertainty  of  T(sun  estimates  is  comparable  for  both  Langley  and  meteorological 
techniques,  even  though  small  differences  are  noted  above  Ka  band. 
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TABLE  II.  Intercomparison  of  Langley  and  meteorological 

ESTIMATES  OF  WEIGHTED  SUN  BRIGHTNESS  TEMPERATURE 


Langley 

Meteorological 

/ 

rGHzl 

Tbsuti[K] 

O T*  [K] 

lBsun L  J 

T-zclr  [Np] 

°V„JNP] 

Tbsuti  [K] 

aT*  [K] 

lBsun L  J 

23.8 

120.82 

0.96 

0.098 

0.005 

119.99 

1.11 

31.4 

182.78 

1.03 

0.043 

0.004 

183.40 

1.47 

72.5 

570.56 

7.19 

0.304 

0.008 

569.07 

7.49 

82.5 

719.22 

10.90 

0.183 

0.010 

704.92 

8.12 

B.  Results  for  path  attenuation  estimate 

The  ST-MWR  technique  is  intrinsically  based  on  a  variable  antenna  pointing  in  order  to  follow  the  Sun 
movement  along  its  ecliptic.  In  order  to  make  atmospheric  path  attenuation  values  intercomparable,  we  can 
show  all  results  in  terms  of  the  zenith-equivalent  path  attenuation,  obtained  by  multiplying  each  estimated 
path  attenuation  for  the  corresponding  air  mass. 

As  an  example,  Fig.  8  (a-d)  shows  the  time  series  of  zenith-equivalent  attenuation  estimates  at  23.8,  31.4, 
72.5  and  82.5  GHz  respectively,  for  the  cloudy  case  of  Fig.  2.  Predictions,  based  on  models  from  PPM  Pol  in 
(26),  PPM  PolDEx  in  (27),  and  Clr  in  (28),  are  compared  against  ST-MWR  attenuation  retrievals.  All 
parametric  prediction  models,  as  expected  in  this  case,  provide  quite  good  results  in  agreement  with  ST- 
MWR  time  series. 

Fig.  9  (a-d)  shows  a  time  series  similar  to  Fig.  8,  but  for  the  intense  rainy  event  of  Fig.3.  Especially  at  72.5 
and  82.5  GHz,  as  expected,  the  application  of  the  Clr  model  is  not  adequate  to  retrieve  path  attenuation 
during  rain. 

For  an  overall  quantitative  estimation  of  the  PPM  prediction  approaches,  Fig.  10  (a-d)  shows  the 
scatterplot  of  PolDEx  and  Clr  model  estimates  for  each  frequency  versus  the  corresponding  ST-MWR  A  for 
all  weather  conditions  and  the  all  available  dataset.  Fig.  11  (4x4)  shows  the  same  comparison  as  in  Fig.  10, 
but  for  Pol  model.  Table  III  quantifies  the  comparison  in  terms  of  RMSE  and  AvE,  where  the  error  is 
defined  as  the  difference  between  prediction  model  under  consideration  and  ST-MWR. 


(a) 

Zenithal  Path  Att.  -  72.5  GHz  -  28-May-2015 


(c) 


Zenithal  Path  Att  -  31.4  GHz  -  28-May-201 5 


(b) 

Zenithal  Path  Att.  -  82.5  GHz  -  28-May-201 5 

8 


<d) 


Fig.  8.  Time  series  of  zenith-equivalent  path  attenuation  estimates  on  May  28,  2015  (cloudy  case)  at  a)  23.8  GHz;  b) 

3 1 .4  GHz;  c)  72.5  GHz;  and  d)  82.5  GHz 
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Zenithal  Path  Att.  -  23.8  GHz  -  30-Jun-2015 


Zenithal  Path  Att.  -  31.4  GHz  -  30-Jun-2015 


(a) 

Zenithal  Path  Att.  -  72.5  GHz  -  30-Jun-2015 


Zenithal  Path  Att.  -  82.5  GHz  -  30-Jun-2015 


(c) 


Fig.  9.  Time  series  of  zenith-equivalent  path  attenuation  estimates  on  June  30,  2015  (rainy  event)  at  a) 
23.8  GHz;  b)  31.4  GHz;  c)  72.5  GHz;  and  d)  82.5  GHz. 


23.8  GHz 


72.5  GHz 


31.4  GHz 


82.5  GHz 


Fig.  10.  Scatterplot  of  ST-MWR  path  attenuation  for  each  frequency  versus  path  attenuation  estimated 
from  PolDEx  and  Clr  models  for  all  weather  condition. 
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23.8  GHz 


31.4  GHz 


72.5  GHz  82.5  GHz 


Fig.  11.  Scatterplot  of  ST-MWR  path  attenuation  for  each  frequency  versus  path  attenuation  estimated 
from  Pol.  for  all  weather  condition.  Results  of  the  Clr  models  are  shown  as  a  reference. 


TABLE  III.  Estimated  zenith-equivalent  attenuation  error 
CONSIDERING  THE  ST-MWR  AVAILABLE  DATASET  IN  20 1 5  IN  ROME,  NY  FOR 
POL  AND  AND  POLDEX  MODELS 


Pol 

PolDex 

Clr-Tmr 

f 

AvE[ 

RMS 

AvE[ 

RMS 

AvE[ 

RMS 

[GHz] 

dB] 

E[dB] 

dB] 

E[dB] 

dB] 

E[dB] 

23.8 

- 

0.250 

- 

0.126 

- 

0.135 

0.0224 

2 

0.0066 

8 

0.0271 

0 

31.4 

- 

0.352 

+0.00 

0.147 

+0.00 

0.219 

0.0068 

6 

18 

5 

81 

7 

72.5 

- 

0.323 

- 

0.383 

- 

0.619 

0.0212 

6 

0.0046 

2 

0.0709 

9 

82.5 

- 

0.428 

- 

0.379 

- 

0.616 

0.1258 

0 

0.0680 

7 

0.1624 

6 

From  Fig.  10  and  11  it  is  worth  noting  that  the  PPM  Pol  prediction  tends  to  underestimate  significantly 
path  attenuation  at  Ka  band,  showing  performances  even  worse  than  Clr  prediction.  On  the  other  hand,  the 
PolDEx  approach  shows  a  fairly  good  correlation  with  reference  ST-MWR  data  for  all  frequencies,  whereas 
the  Clr  model  is  much  worse  than  the  other  models  at  V  and  W  band.  Note  that  the  latter  is  not  even 
applicable  in  rainy  conditions  and  at  higher  frequencies  when  the  atmosphere  is  opaque.  Table  III  confirms 
these  visual  considerations  stressing  the  very  low  bias  of  PolDEx  estimates  and  its  flexibility  in  dealing  with 
both  clear  and  rainy  cases.  Note  that  the  scores  of  the  Clr  model  are  slightly  optimistic  as  they  computed  for 
a  reduced  number  of  applicable  points. 

In  (30),  coupled  with  (27),  we  have  introduced  a  dual-frequency  approach  to  path  estimation  retrieval. 
Table  IV,  similarly  to  Table  III,  quantifies  the  results  of  this  reduced-channel  techique  in  terms  of  RMSE 
and  AvE.  Note  how  path  attenuation  retrieval  errors  from  Pol2DEx  model  is  still  relatively  small; 
nevertheless,  higher  accuracy  is  obtained  using  more  than  2  Ka-band  channels. 


TABLE  IV.  Same  as  Table  III,  but  for  the  dual-frequency  prediction 

MODEL 


Pol  (bifreq) 

PolDex  (bifreq) 

Tmr 

/ 

AvE[d 

RMSE 

AvE[d 

RMSE 

AvE[d 

RMSE 

[GHz] 

B] 

[dB] 

B] 

[dB] 

B] 

[dB] 

23.8 

- 

0.242 

- 

0.127 

- 

0.135 

0.0178 

5 

0.0026 

2 

0.0271 

0 

31.4 

0.004 

0.351 

0.013 

0.147 

+0.00 

0.219 

5 

6 

4 

1 

81 

7 
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C.  Error  budget  analysis 

Sect.  II.C  has  been  devoted  to  the  discussion  of  the  theoretical  error  analysis  of  path  attenuation  estimates 
with  respect  to  several  error  sources.  By  using  previous  results,  obtained  from  the  AFRL-MWR 
measurements,  we  can  now  evaluate  in  a  more  quantitative  way  the  overall  error  budget.  The  latter  is  shown 
in  Table  V  which  summarizes  the  expected  errors  in  A  and  at  23.8,  31.4,  72.5  and  82.5  GHz  with  respect  to 
the  uncertainty  of  the  various  sources. 


TABLE  V.  ERROR  EVALUATION  FOR  PATH  ATTENUATION  ESTIMATES  USING 
AFRL-MWR  MEASUREMENTS 


8 A  with  respect  to  8  A TA 

SA  with  respect  to  8Tgsun 

/[GHz]  j 

SArqK 

1 

54  [dB] 

/[GHz] 

srBsun[ 

K] 

54  [dB] 

23.8 

4 

-0.20 

23.8 

0.83 

0.030 

31.4 

5 

-0.13 

31.4 

0.62 

0.015 

72.5 

12 

-0.13 

72.5 

1.49 

0.011 

82.5  1 

19 

-0.14 

82.5 

14.3 

0.085 

Sensitivity  SA  with  respect  to  SATa  is  given  by  (22a);  on  the  left  side  of  Table  V  reports  the  attenuation 
uncertainties  for  typical  SATa  in  clear  air.  The  latter  are  derived  from  the  variability  of  AFRL-MWR  data 
during  the  ooS  and  twS  switching.  In  this  case  5A  are  less  than  0.2  dB  at  all  frequencies,  slightly  higher  at  Ka 
band.  On  the  right  side  of  Table  V,  a  the  sensitivity  analysis  of  SA  is  shown  with  respect  to  uncertainties 
STzsun-  The  latter  has  been  derived  from  the  differences  between  Langley  and  meteorological  estimates 
shown  in  Table  V.  Errors  in  path  attenuation  retrievals  are  less  than  0.1  dB  at  all  frequencies.  Finally,  we 
have  estimated  the  error  in  assuming  the  horizontal  homogeneity  in  clear  sky  through  the  analysis  of  the  ooS 
time  series  at  the  same  elevation  and  different  azimuths.  Attenuation  uncertainties  are  estimated  to  be  0.017, 
0.007,  0.023,  and  0.027  dB  for  an  average  azimuth  distance  of  5°.  As  such,  the  horizontal  homogeneity 
assumption  holds. 

V.  Conclusion 

Sun-tracking  microwave  radiometry  has  been  introduced  to  estimate  the  atmospheric  path  attenuation  in 
all-weather  conditions  at  Ka,  V  and  W  band.  A  detailed  theoretical  framework  has  been  proposed  to  describe 
the  Sun-tracking  microwave  radiometric  measurement  modes  and  to  evaluate  the  overall  error  budget  with 
respect  to  several  sources  of  uncertainties.  This  approach  has  clearly  identified  the  critical  assumptions 
behind  the  ST-MWR  data  processing  such  as  the  accurate  knowledge  of  the  antenna  beam  width,  the 
uniformity  of  incident  brightness  temperature  within  the  antenna  beamwidth  and  the  atmospheric  stationarity 
within  each  ST  system  switch. 

The  weighted  brightness  temperature  of  the  Sun  Tgsun  has  been  estimated  by  means  of  two  different 
techniques,  based  on  the  elevation-scanning  Langley  method  and  the  surface  meteorological  data  method. 
These  two  techniques  show  comparable  results,  but  the  first  one  needs  a  careful  selection  of  candidate  clear- 
air  days  whereas  the  second  one  is  depending  on  the  external  weather  station  data  and  daily  variability  of 
clear  air  extinction.  The  ST-MWR  methodology  has  been  applied  to  AFRL  available  radiometric 
measurements  at  Ka,  V  and  W  band  in  Rome  (NY,  USA)  during  2015,  in  order  to  test  two  new  physically- 
based  prediction  models  for  path  attenuation  estimation.  The  single-frequency  double-exponential  parametric 
model  seems  to  outperform  the  multi-frequency  polynomial  model  when  compared  to  ST-MWR  retrievals. 
The  clear-air  based  approach  to  estimate  path  attenuation,  as  expected,  is  unsuitable  to  estimate  path 
attenuation  in  heavy  cloudy  and  rainy  conditions  and  at  higher  millimeter-wave  frequency. 

These  results  show  an  appealing  potential  of  ST-MWR  technique,  which  can  be  exploited  to  overcome  the 
overall  costs  and  the  logistic  difficulty  to  accomplish  satellite-to-Earth  radiopropagation  experiments  in  the 
unexplored  millimeter- wave  and  submillimeter-wave  frequency  region.  A  further  validation  of  the  ST-MWR 
technique  should  foresee  by  enlarging  the  available  dataset  and  by  comparing  with  collocated  measurements 
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performed  by  a  satellite  beacon  receiver,  even  though  the  telecommunication  bandwidths  are  typically 
different  from  the  available  MWR  ones.  The  deployment  of  several  AFRL  MWRs,  where  PPM  are  applied  to 
estimate  path  attenuation  after  a  verification  through  ST-MWR,  might  allow  very  promising  site  diversity 
experiments  at  V  and  W  band. 


Appendix.  Parametric  prediction  models 

This  appendix  provides  the  tables  of  the  regression  coefficients  and  modeling  functions  used  to  express  the 
physically-based  parametric  models  (PPMs),  discussed  in  sect.  III.B.  Two  vPPM  ersions  are  considered:  i) 
multifrequency,  assuming  to  have  a  disposal  all  4  AFRL-MWR  channels  (sect.  A);  ii)  dual-frequency, 
assuming  to  have  only  Ka  band  channels  at  23.8  and  31.4  GHz  (sect.  B). 

A.  Multifrequency  PPM  model 

Table  A1  provides  the  regression  coefficients  expressions  to  apply  Pol  model  in  (26),  whereas  Table  A2  the 
same  but  for  the  DEx  model  in  (27b)  for  each  considered  frequency / at  23.8,  31.4,  72.5  and  82.5  GHz.  Table 
A3  provides  the  Pol  and  PolDEx  model  intercomparison  in  terms  of  error  indexes. 


TABLE  A 1 .  PPM  POLYNOMIAL  MODEL  COEFFICIENTS  IN  (26) 

23.8  GHz  ~ 

a,  =  +2.51  ■  10"2  -  1.42  ■  1CT2  ■  m  +  2.31  ■  10“3  ■  m2 

a2  =  +8.16  ■  10“3  -  4.12  ■  10“3  ■  m  +  5.03  ■  10“s  ■  m2 

a3  =  -4.18  ■  10-3  +  2.97  ■  10“3  ■  m  -  3.61  ■  10"4  ■  m2 

a4  =  +1.84  ■  10“3  -  1.65  ■  10“3  ■  m  +  5.81  ■  10"4  ■  m2 

b4  =  +6.04  ■  10"s  -  3.19  ■  10“s  ■  m  +  6.32  ■  10“6  ■  m2 

b2  =  -2.04  ■  10“s  +  2.46  ■  10“s  ■  m  -  3.07  ■  10“6  ■  m2 

b3  =  +2.67  ■  10“s  -  2.34  ■  10“s  ■  m  +  3.60  ■  10“6  ■  m2 

b4  =  -1.95  ■  10“s  +  1.60  ■  10“s  ■  m  -  3.04  ■  10“6  ■  m2 

31.4  GHz  ~ 

a,  =  -3.37  ■  10"3  +  1.99  ■  10“3  ■  m  -  4.59  ■  10"4  ■  m2 

a2  =  +4.04  ■  10-2  -  2.03  ■  10“2  ■  m  +  2.03  ■  10“3  ■  m 2 

a3  =  -6.13  ■  10“3  +  4.07  ■  10“3  ■  m  -  4.27  ■  10“s  ■  m2 

a4  =  +9.42  ■  10“4  -  1.03  ■  10“3  ■  m  +  6.19  ■  10"4  ■  m2 

b4  =  +1.60-  10“s  -  1.71-  10“s  -m  + 4.11-  10“6  ■  m2 
b2  =  +2.84  ■  10“s  +  1.24  ■  10“s  ■  m  +  1.55  ■  10“7  ■  m2 

b3  =  +5.53  ■  10“s  -  5.71  ■  10“s  ■  m  +  7.54  ■  10“6  ■  m2 

b4  =  -3.81  ■  10“s  +  3.62  ■  10“s  ■  m  -  5.94  ■  10“6  ■  m2 

72.5  GHz  ~ 

a4  =  +3.37  ■  10"3  -  5.86  ■  10“3  ■  m  +  9.66  ■  10"4  ■  m2 

a2  =  +2.94  ■  10“2  -  4.11  ■  10“2  ■  m  +  7.13  ■  10“3  ■  m2 

a3  =  +2.75  ■  10“2  -  1.91  ■  10“2  ■  m  +  3.26  ■  10“3  ■  m2 

a4  =  -1.76  ■  10“2  +  2.01  ■  10“2  ■  m  -  3.33  ■  10“3  ■  m2 

b4  =  -7.61-  10“s  +  6.42-  10“s  -m-  1.11-  10“s  ■  m2 
b2  =  +4.00  ■  10"4  -  8.52  ■  10-s  ■  m  +  1.48  ■  10“s  ■  m 2 

b3  =  +2.51  ■  10“4  -  1.29  ■  10"4  ■  m  +  2.25  ■  10“s  ■  m2 

b4  =  -1.17  ■  10"4  +  4.11  ■  10“s  ■  m  -  7.44  ■  10“6  ■  m2 

82.5  GHz  ~ 

aj  =  +9.71  ■  10"3  -  7.15  ■  10-3  ■  m  +  1.18  ■  10“3  ■  m2 

a2  =  +4.10  ■  10“2  -  5.28  ■  10“2  ■  m  +  9.06  ■  10“3  ■  m2 

a3  =  -1.19  ■  10-2  +  9.29  ■  10“3  ■  m  -  1.65  ■  10“3  ■  m2 

a4  =  +1.40  ■  10“2  -  2.71  ■  10“3  ■  m  +  6.23  ■  10"4  ■  m2 

b4  =  -1.60  ■  10"4  +  1.13  ■  10"4  ■  m  -  1.95  ■  10“s  ■  m2 

b2  =  +4.16  ■  10-4  -  7.08  ■  10-5  ■  m  +  1.23  ■  10“s  ■  m2 

b3  =  +2.05  ■  10"4  -  1.40  ■  10"4  ■  m  +  2.47  ■  10“s  ■  m2 

b4  =  -4.81  ■  10“s  +  3.25  ■  10“s  ■  m  -  6.09  ■  10“6  ■  m2 
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TABLE  A3.  PPM  MODEL  INTERCOMPARISON 


23.8  GHz 


Pol 

PolDex 

Elev 

AvE[d 

RMSE 

IA 

AvE[d 

RMSE 

IA 

90° 

-0.213 

1.495 

0.801 

-0.050 

0.507 

0.986 

75° 

-0.216 

1.557 

0.798 

-0.047 

0.511 

0.987 

On 

O' 

-0.229 

1.777 

0.784 

-0.042 

0.558 

0.988 

45° 

-0.283 

2.293 

0.748 

-0.072 

0.838 

0.981 

36° 

-0.386 

2.906 

0.705 

-0.176 

1.398 

0.958 

30° 

-0.527 

3.569 

0.662 

-0.338 

2.144 

0.918 

25° 

-0.698 

4.380 

0.620 

-0.534 

2.992 

0.874 

20° 

-0.731 

5.509 

0.596 

-0.546 

3.623 

0.881 

31.4  GHz 

Pol 

PolDex 

Elev 

AvE[d 

RMSE 

IA 

AvE[d 

RMSE 

IA 

90° 

-0.317 

2.244 

0.608 

-0.156 

1.069 

0.950 

75° 

-0.324 

2.336 

0.601 

-0.153 

1.100 

0.950 

60° 

-0.354 

2.657 

0.576 

-0.150 

1.222 

0.951 

45° 

-0.443 

3.394 

0.511 

-0.167 

1.556 

0.946 

36° 

-0.580 

4.257 

0.435 

-0.223 

1.965 

0.940 

30° 

-0.757 

5.182 

0.361 

-0.300 

2.344 

0.938 

25° 

-0.976 

6.318 

0.288 

-0.376 

2.708 

0.943 

20° 

-1.096 

7.940 

0.237 

-0.561 

4.842 

0.845 

72.5  GHz 

Pol 

PolDex 

Elev 

AvE[d 

RMSE 

IA 

AvE[d 

RMSE 

IA 

90° 

-0.391 

2.892 

0.834 

-0.468 

2.694 

0.863 

75° 

-0.392 

3.037 

0.828 

-0.454 

2.798 

0.861 

60° 

-0.416 

3.554 

0.803 

-0.421 

3.164 

0.856 

45° 

-0.582 

4.802 

0.739 

-0.441 

4.020 

0.843 

36° 

-0.953 

6.331 

0.660 

-0.623 

5.062 

0.825 

30° 

-1.494 

8.018 

0.580 

-0.957 

6.252 

0.805 

25° 

-2.176 

10.084 

0.503 

-1.405 

7.807 

0.778 

20° 

-2.383 

12.780 

0.454 

-1.383 

10.154 

0.742 

82.5  GHz 

Pol 

PolDex 

Elev 

AvE[d 

RMSE 

IA 

AvE[d 

RMSE 

IA 

90° 

-0.404 

2.988 

0.829 

-0.458 

2.903 

0.838 

75° 

-0.407 

3.135 

0.822 

-0.451 

3.021 

0.836 

60° 

-0.436 

3.655 

0.799 

-0.437 

3.441 

0.826 

45° 

-0.602 

4.901 

0.739 

-0.490 

4.438 

0.801 

36° 

-0.958 

6.417 

0.664 

-0.688 

5.641 

0.771 

30° 

-1.481 

8.083 

0.590 

-1.013 

6.959 

0.742 

25° 

-2.162 

10.121 

0.519 

-1.430 

8.591 

0.714 

20° 

-2.501 

12.813 

0.482 

-1.419 

10.860 

0.700 
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TABLE  A2.  PPM  DOUBLE  EXPONENTIAL  MODEL  COEFFICIENTS  IN  (21  A) 

23.8  GHz 

a  =  +1.26  -  9.27  ■  10"1  ■  m  +  1.78  ■  10"1  ■  m2 
b  =  +6.85  ■  10“3  +  7.61  ■  1CT3  ■  m  -  2.00  ■  1CT3  ■  m2 
c  =  +2.18  ■  10"14  +  2.92  ■  lO"30  ■  m  -  2.70  ■  1(T30  ■  m2 
d  =  +1.17  ■  10"1  +  0.00  ■  m  -  1.42  ■  lO"17  ■  m2 
h  =  -2.01  ■  10"1  +  6.18  ■  10"1  ■  m  -  1.44  ■  10"1  ■  m 2 

31.4  GHz 


a  =  +1.01  -  5.38  ■  10"1  ■  m  +  7.38  ■  lO"2  ■  m2 
b  =  +1.57  ■  10"2  -  4.56  ■  10“3  ■  m  +  1.83  ■  10"3  ■  m2 
c  =  -9.36  ■  10"14  +  1.57  ■  lO"13  ■  m  -  4.17  ■  10_14  ■  m2 
d  =  +1.25  ■  10"1  -  1.12  ■  10“2  ■  m  +  3.14  ■  1(T3  ■  m2 
h  =  +5.34  ■  10"1  -  1.46  ■  10"1  ■  m  +  3.26  ■  lO"2  ■  m2 


72.5  GHz 


a  =  +1.75  -  1.11  ■  m  +  2.03  ■  10"1  ■  m2 
b  =  +7.56  ■  10"3  +  2.46  ■  10“3  ■  m  -  4.31  ■  10"4  ■  m2 
c  =  -1.85  ■  10“9  -  3.88  ■  10"2S  ■  m  -  7.07  ■  10"26  ■  m2 
d  =  +8.00  ■  10“2  -  1.30  ■  10"17  ■  m  -  9.49  ■  10"18  ■  m2 
h  =  +1.42  -  1.08  ■  m  +  2.37  ■  10"1  ■  m2 


82.5  GHz 


a  =  +1.62  -  9.43  ■  10"1  ■  m  +  1.68  ■  10"1  ■  m2 
b  =  +8.85  ■  10"3  +  1.06  ■  10“3  ■  m  -  1.60  ■  10"4  ■  m2 
c  =  +3.66  ■  lO"10  -  2.95  ■  10"12  ■  m  +  4.76  ■  10“13  ■  m2 
d  =  +8.51  ■  10"2  -  1.30  ■  10"17  ■  m  -  4.75  ■  10"18  ■  m 2 
h  =  +9.06  ■  10"1  -  6.27  ■  10"1  ■  m  +  1.25  ■  10"1  ■  m2 


B.  Dual-frequency  PPM 

Table  A4  provides  the  regression  coefficients  expressions  to  apply  Pol2  model  in  (30). 


TABLE  A4.  DUAL-FREQUENCY  POLYNOMIAL  MODEL  COEFFICIENTS 
23.8  GHz 

Qj  =  +2.02  ■  10"3  -  1.14  ■  10"2  ■  m  +  2.02  ■  10"3  ■  m2 
a2  =  +3.86  ■  10“3  -  2.82  ■  10“3  ■  m  +  3.09  ■  10"4  ■  m2 
=  +9.62  ■  10"s  -  6.71  ■  10"s  ■  m  +  1.28  ■  10"s  ■  m2 
b2  =  -1.04  ■  10"s  +  2.12  ■  10"s  ■  m  -  3.65  ■  10"6  ■  m2 
31.4  GHz 

at  =  -6.19  ■  10"3  +  5.84  ■  10"3  ■  m  -  8.77  ■  10"4  ■  m2 
a2  =  +3.19  ■  10"2  -  2.02  ■  10"2  ■  m  +  3.23  ■  10“3  ■  m2 
=  +3.45  ■  10"s  -  4.32  ■  10"s  ■  m  +  8.55  ■  10"6  ■  m2 
b2  =  +7.25  ■  10"s  -  1.26  ■  10"s  ■  m  +  2.36  ■  10"6  ■  m2 
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